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Abstract. Eternal inflation arising from a potential landscape predicts that our universe is 
one realization of many possible cosmological histories. One way to access different cosmo¬ 
logical histories is via the nucleation of bubble universes from a metastable false vacuum. 
Another way to sample different cosmological histories is via classical transitions, the cre¬ 
ation of pocket universes through the collision between bubbles. Using relativistic numeri¬ 
cal simulations, we examine the possibility of observationally determining if our observable 
universe resulted from a classical transition. We find that classical transitions produce spa¬ 
tially infinite, approximately open Friedman-Robertson-Walker universes. The leading set 
of observables in the aftermath of a classical transition are negative spatial curvature and 
a contribution to the Cosmic Microwave Background temperature quadrupole. The level of 
curvature and magnitude of the quadrupole are dependent on the position of the observer, 
and we determine the possible range of observables for two classes of single-scalar field mod¬ 
els. For the first class, where the inflationary phase has a lower energy than the vacuum 
preceding the classical transition, the magnitude of the observed quadrupole generally falls 
to zero with distance from the collision while the spatial curvature grows to a constant. For 
the second class, where the inflationary phase has a higher energy than the vacuum pre¬ 
ceding the classical transition, the magnitude of the observed quadrupole generically falls to 
zero with distance from the collision while the spatial curvature grows without bound. We 
find that the magnitude of the quadrupole and curvature grow with increasing centre of mass 
energy of the collision, and explore variations of the parameters in the scalar field lagrangian. 
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1 Introduction 

In the picture of eternal inflation (see Refs. [1-3] for a review), our observable universe is 
part of a multiverse in which different regions of space are occupied by different field values 
in a potential landscape (for example, the simple potential landscape shown in Eig. 1). 
In eternal inflation, a region of the universe initially occupying a false vacuum gives rise to 
bubble universes that ultimately reach different vacua after a period of cosmological evolution. 
If the rate of bubble formation does not outpace the background expansion rate, inflation 
becomes eternal and the process of bubble formation carries on forever. 

Although the phase transition out of the false vacuum does not complete, all of the bub¬ 
bles that form during eternal inflation undergo collisions with others [4, 5]. If our Observable 
Universe is contained within one such bubble, collisions provide a possible observational test 
of eternal inflation [6]. A substantial body of work has been devoted to determining the 
outcome of bubble collisions and predicting the observational signatures [5-46] , and the first 
searches for the predicted pattern in the cosmic microwave background (CMB) anisotropy 
have put constraints on the number and amplitude of collision signatures in our observable 
Universe [10-12, 43]. 

In a simple potential landscape with only two vacua, colliding bubbles merge. When the 
scalar field potential has three or more local minima, bubble collisions can produce domain 
walls separating regions occupying different vacua. In some cases, the collision can produce a 
region that occupies a new vacuum, distinct from the vacua in either of the colliding bubbles. 
Such an event is known as a classical transition [32], and is an alternative to quantum me¬ 
chanical bubble nucleation as a mechanism for populating vacua and cosmological histories in 
a potential landscape. In some cases, classical transitions may represent the only mechanism 
for accessing certain regions of a potential landscape [13]. If our Universe inhabits such a 
corner of the potential landscape, then it must have arisen from a classical transition. In 
this case, the rate at which classical transitions occur is irrelevant, since bubble collisions are 
guaranteed during eternal inflation. However, in the absence of a more complete theory, it 
is difficult to judge how generic classical transitions might be. 
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As we will see below, the Universe to the future of a classical transition is approximately 
Friedman-Robertson-Walker (FRW) with negative curvature. In analogy to models of open 
inflation [47], one can invoke a period of slow-roll inflation after the classical transition occurs 
in order to dilute this spatial curvature and produce the appropriate density perturbations 
to seed large scale structure. In this way, one can create models where our Universe might 
exist to the future of a classical transition. The aim of this paper is to investigate the 
possibility that our Observable Universe could have been produced by a classical transition, 
and compute the observational signatures of this scenario. 

Because classical transitions are an intrinsically non-linear phenomenon, it is necessary 
to employ numerical simulations. To determine the cosmology to the future of a classical 
transition, we must perform these simulations in full General Relativity. Fortunately, the 
collision spacetime possesses a hyperbolic symmetry which allows us to simulate the entire 
collision spacetime using a simulation in one space and one time dimension. This symmetry 
makes performing simulations computationally cheap, allowing us to simulate many different 
cases. Throughout the paper, we use the simulation framework developed in Refs. [16, 17, 48]. 
This framework directly links the scalar field landscape underlying eternal inflation to the 
predictions for cosmological observables, allowing for observational constraints to be placed 
on the theory underlying eternal inflation. 

We find that classical transitions produce spatially infinite approximately FRW uni¬ 
verses. The leading cosmological observables are negative spatial curvature and an approxi¬ 
mately quadratic comoving curvature perturbation with planar symmetry. Such a curvature 
perturbation maps to a temperature quadrupole in the Cosmic Microwave Background radi¬ 
ation. The prediction for curvature and the CMB quadrupole are correlated, and depend on 
the distance of the observer (in an appropriate set of coordinates) from the spatial location of 
the collision that caused the classical transition. In addition to observer position, the initial 
separations between the two colliding bubbles is a random variable. In addition to depending 
on the scalar field lagrangian, the observables also depend on the initial separation. There¬ 
fore, any specific scalar field lagrangian giving rise to classical transitions yields an ensemble 
of predictions for observation determined by the initial separation of the colliding bubbles 
and the position of an observer. We classify the range of predictions allowed in the context 
of two classes of models for the scalar field potential. 

The plan of the paper is as follows. In Sec. 2 we review the mechanism of classical 
transitions and the expected geometries allowed in GR. In Sec. 3 we describe the simulation 
framework and methods for extracting cosmological observables. In Sec. 4 we present the 
results from our simulations, and in Sec. 5 we conclude 

2 Classical Transitions 

In Fig. 1, we outline a number of properties of classical transitions. Consider a potential with 
local minima A, B, and C. We will consider cases where the potential Vc at C is both lower 
(dashed) and higher (solid) than the potential Ug at S. Vacuum A is a false vacuum and 
bubbles of B spontaneously form through quantum tunnelling, mediated by the Coleman and 
De Luccia (CDL) [49-51]. Many collisions between bubbles of B occur, in spite of the fact 
that the phase transition does not complete in eternal inflation. For pairwise collisions, which 
are most relevant in the regime of small nucleation rates, we are able to work in the centre 
of mass frame, where the two colliding bubbles nucleate at the same coordinate time [7]. In 
the absence of gravity, the lorentz factor of the bubble walls at the time of collision is given 
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by: 


( 2 . 1 ) 


Ax 

where R is the initial radius of the colliding bubbles and Ax is the proper distance between 
the bubbles in the centre of mass frame. The lorentz factor characterizes the centre of mass 
energy of the collision. 

When the bubbles of B collide, the field in the immediate aftermath of the collision is 
displaced. If the centre of mass energy is large enough (7 ^ 0(1)), this displacement reaches 
a maximum size, given by the superposition of two bubble profiles [33]. In the limit of 
infinite boost, superposition holds exactly; this is known as the free passage approximation. 
For vacuum bubbles, as we consider here, the amplitude of each of the bubble profiles, and 
therefore the maximum field excursion, is the width of the potential barrier, Ac/), separating 
A and B. When the structure of the potential is such that the field is pushed into a local 
minimum, as in Fig. 1, a classical transition results, producing a spacetime region containing 
phase C. 

The ultimate fate of the spacetime region containing C depends on the vacuum energies 
of the three phases. This is depicted in Fig. 1 . For Vb > Vc, region of C expands due to 
the outward pressure gradient across the domain wall separating C from B. The result is a 
lasting region of (7; this is known as a Normal geometry. Without gravity, when Vb < Vc, ^.n 
Oscillatory geometry is produced. In an Oscillatory geometry, an initially expanding pocket 
of C re-collapses due to the inward pressure gradient across the domain wall separating C 
from B. The collision makes a new region of C which itself re-collapses, and this process 
repeats until the colliding walls dissipate into scalar radiation. No lasting region of C is 
produced in an Oscillatory geometry. Including gravitational effects [14, 52], it is possible 
to produce a lasting region of C even when Vb < Vc- These geometries are referred to 
as Marginally Repulsive, and occur when the gravitationally repulsive effects of the domain 
walls separating C from B are strong enough to prevent the region of C from re-collapsing. 
This can occur when Vb with Vc bounded by the constraint ^ 


Vc< 


(^4 + Vb? 

4Vk 


( 2 . 2 ) 


In this paper, we focus on Normal and Marginally Repulsive geometries, both of which 
contain lasting regions of C. 

Previous work has characterized and simulated classical transitions [13, 14, 30-33]. The 
focus of the present work is on determining what the cosmological observables would be if 
we inhabited a region produced by a classical transition, a question to which we now turn. 


3 Simulation Framework 


3.1 Scalar field potential 

We use a quadratic inflationary potential with two Gaussian bumps for our model potential 
landscape: 


V((/)) = AiExp 


2A(/)f_ 


± A2EXP 


(T 


2 A (/)2 


'2 J 


+ -m 


- 


(3.1) 


^Note that Eq. 2.2 implies that Vc < Va for a marginally repulsive geometry. Therefore, one can never 
attain a potential higher than the False vacuum by a classical transition. 
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Vb > Vc Vb « Vc Vb < Vc 



Figure 1. A scalar field potential that gives rise to classical transitions (top) and the possible 
classical transition geometries (bottom). The scalar field potential has three local minima: A, 
and C. We will consider cases where the potential Vc at C is both lower (dashed) and higher (solid) 
than the potential Vb at B. The field excursion between vacuum A and vacuum B is A0. In the free 
passage approximation, the distance the field is pushed in the future of a collision between two bubbles 
of B embedded in a background of A is bounded by Acj). Including gravity, there are a number of 
distinct classical transition geometries (bottom). We consider bubbles of B nucleated in a background 
false vacuum of A. Bubbles of B form with an initial radius i?, and in the centre of mass frame, are 
separated by a proper distance Ax. For a Normal geometry, which corresponds to Vb > Vc, the 
region of C produced by the collision expands indefinitely. In an Oscillatory geometry, corresponding 
to Vc Vb, the initially expanding region of C eventually collapses, spawning a new region of C 
which expands and contracts. Eventually this oscillation terminates in the release of scalar radiation, 
and vacuum C is no longer present in the spacetime. Repulsive and Marginally Repulsive geometries 
occur when Vb and Vc are comparable. For Vb ^ Vc, 8i Marginally Repulsive geometry results, in 
which the region of C expands but is prevented from contracting by the gravitational field of the 
domain wall and the positive energy of vacuum C. This gives rise to a lasting region of C. 


In Eq. 3.1, “ + ” indicates whether the second barrier is above (AIP) or below (BIP) the 

inflationary plateau respectively. The inflationary component of the potential is necessary to 
dilute the curvature of spatial slices in the classical transition spacetime, as we will see below. 
From the discussion in the previous section, the AIP models produce classical transitions 
with Normal geometries and the BIP models produce Marginally Repulsive geometries. This 
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0/Mp 


Figure 2. The potential Eq. 3.1 for a characteristic choice of parameters. The true vaccum is located 
at 0 = 2Mp, the solid line is the AIP model; the dashed line is the BIP model. 


choice of potential has no microphysical motivation, but has enough freedom to explore the 
phenomenology of classical transitions. A(/)i /2 the widths of the bumps, a is the shift 
of the second bump from the first one, and 00 specifies the location of the true vacuum. 
Following Ref. [17], we choose to parameterize Ai and A 2 by: 

Ai = /3i|m^0oA0ie^/^|, A 2 = /32|^^0oA02e^^^| (3.2) 

In general, there are five free variables associated with the potential Eq. 3.1: the barrier 
widths (A0 x/ 2)5 the parametrized bump heights {P 1 / 2 ) ^^nd the shift (a). Below, we will 
measure A 0 x /2 ^ terms of Mp] /3x/2 is dimensionless. 

We first discuss the parameters of the AIP model. We can make a few simplifications 
to reduce the dimensionality of parameter space. The free passage approximation suggests 
that the structure of the second barrier (characterized by A02 sind 02 ) is largely irrelevant in 
the AIP model. Free passage also suggests that there will be large parameter degeneracies, 
as the most important determinants of the outcome of a collision are the initial profile of 
the bubble wall (given by the analytic continuation of the CDL instanton) and the initial 
separation between bubbles. 

In Fig. 3, we show the potential for the AIP model for varying 0i, A0i, and cr; the 
associated instanton profiles are shown in Fig. 4. We hold 02 and A02 fixed for the AIP 
model. However, note that the first barrier can affect the shape of the second barrier in some 
regions of parameter space. 

In light of the free passage approximation, we can predict that the degeneracies in 
parameter space are associated with parameter combinations that give nearly identical in¬ 
stanton profiles. As a crude model, we can approximate instanton profiles by their amplitude 
0 amp 5 defined as the distance in field space between the instanton endpoints, and their size i?, 
defined as the position where the field attains the midpoint between the instanton endpoints. 
In the analysis below, we choose combinations of the parameters 0i, A0i, and a that give 
instanton profiles with a constant radius but varying amplitude, and vice versa. We can pair 
A01, 01 and a in three different combinations (A0i — 0i, A0i — a and 0i — a) to vary the 
radius and amplitude. However, the window we are able to simulate in the A0i — 0i sector 
is too narrow to show significant differences in i?, and so we omit this sector in the analysis 
below. 
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Figure 3. Potential Eq. 2 with varying parameters Pi (left panel), Api (centre panel), and a (right 
panel). 
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Figure 4. Instanton profiles for the potentials shown in Fig 3: varying Pi (left panel), Api (centre 
panel), and a (right panel). 


For the BIP model, we are primarily interested in exploring the observables in Marginally 
Repulsive geometries. To this end, from Eq. 2.2, we wish to change the relative vacuum 
energies in such a way to produce a Marginally Repulsive geometry. Varying only /32, as 
shown in Fig. 5, there will be a very small range of values for P 2 that produce a marginally 
repulsive geometry (recall that we must satisfy the inequality Eq. 2.2). Simultaneously 
varying a — P 2 oi A(j )2 — P 2 can produce a variety of cases in which a marginally repulsive 
geometry is produced, since varying a and/or A(j)2 will shift the energy density at the onset 
of slow-roll. 

3.2 Collision Spacetime 

Individual CDL bubbles possess an 50(3,1) symmetry [49-51]. The interior of each bubble 
is described by an open FRW metric, and its wall expands asymptotically to the speed of 
light. We will consider pair-wise collisions between bubbles; in the eternally inflating regime, 
this is sufficient to describe the region of spacetime causally accessible to most observers. 
The spacetime describing the collision between two CDL bubbles has an 50(2,1) symmetry, 
and the metric can, in full generality, be written as [16]: 

Hpds^ = x)dN‘^ + a^{N^ x) cosh^ Ndx‘^ + sinh^ N{dx^ + sinh^ x d^‘^) (3.3) 

where Hp is the false vacuum Hubble parameter, N measures the numbers of e-foldings in 
the false vacuum, x is the direction between the two bubble centers, and x? '0 parameterize 
the hyperbolic symmetry of the collision spacetime. 
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Figure 5. for the BIP model at increasing values of (32^ from P 2 = ^ (blue) to /32 = 9.87 

(yellow). The depth of the intermediate potential minimum increases with increasing /32. 


Due to the hyperbolic symmetry, the collision spacetime can be obtained by doing simu¬ 
lations in one time and one space direction. The initial conditions are specified by the analytic 
continuation of the CDL instanton, and the unknown metric functions a{N^ x) and a{N^ x) as 
well as the scalar field (/)(A^, x) evolved using the Einstein and scalar field equations. Impor¬ 
tantly, the only free parameter in the initial conditions is the initial separation between the 
bubbles in the centre of mass frame. A complete description of the procedure for determining 
the initial value problem and subsequent evolution is described in Refs. [14, 16, 17]. 

3.3 Extracting CMB Observables from the Simulation 

The metric Eq. 3.3 and field 0(A^, x) fully specify the collision spacetime. To make connection 
with cosmological observables, it is necessary to go from this global description to local 
descriptions of the Universe that would be seen by an ensemble of observers. There are a 
variety of approaches to this problem [16, 17, 48, 53]. We use the method outlined in Ref. [48], 
which we briefly summarize here. 

• In the comoving gauge, the scalar field is homogeneous. Therefore, we can extract 
spatial slices in the comoving gauge by finding the surfaces of constant field. We 
choose a reference field value of 00 = 0.15 Mp, which is sufficiently far into inflation 
that the comoving curvature perturbation is frozen in on large scales. 

• To coordinatize the surface of constant field, we use the proper distance u from a 
fiducial point along the surface of constant field, and retain the hyperbolic symmetry 
of the collision spacetime. Projecting the metric Eq. 3.3 yields the spatial metric: 

Hpds\ — du^ + sinh^ N{u){dx^ + sinh^ (3.4) 

• We then define a reference point xq, and make the linear coordinate transformation 

u-u^^ ao(C - Co) (3.5) 

which takes us to the perturbed hyperbolic slicing of FRW: 

Hpds^ — aQ[d^‘^ + (1 ~ 2S) cosh^ + sinh^ pd(p‘^)] (3.6) 
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We can identify (^o as the coordiante associated with xq and ao as the scale factor at 
the reference point. 


We determine the scalar curvature of the metric Eq. 3.6, and integrate the definition 
of the comoving curvature perturbation: 

= -6 + 4V27^(^|xo). (3.7) 

with the boundary conditions 

7^(eoko)=0, ae7^(eoko) = 0. (3.8) 


Neither TZ{^q\xq) nor are physical observables, and so we are free to set 

them to zero. 

We then transform to the anisotropic hyperbolic slicing of FRW introduced in Ref. [16] 
Hpds^ = aQ(l — 27^)+ cosh^ iidp^ + sinh^ pd(p‘^)] (3.9) 


In the vicinity of the reference point we can make the approximation that C — ^5 
the leading correction is given by 




5|77(^o|a;o) 

6 




(3.10) 


which will be small in the window corresponding to our Observable Universe where 
C ~ Co 1- In the following, we will substitute C ^ C- 


In the anisotropic hyperbolic coordinates, we perform a simple translation taking the 
point Co to the origin of coordinates C = 0- 


Finally, we can find the comoving curvature perturbation as a function of Cartesian 
coordinates 

Hlds^ = al ~ ^ [dX‘^ + dY^ + dZ^] (3.11) 

0 “ 4R2 ) 

‘"curv 

where R = In the limit of small C, along the C-direction we have 

i?^ = i?Lv(C^ + 0(C")) (3.12) 

and 

X ~ i?curvCcOsd (3.13) 

We therefore identify C = 1 with the curvature radius Rcurv The curvature radius can 
be written in terms of the present energy density in curvature, as: 

^curv = ~ ~r= (3.14) 

where ido is the present day Hubble constant. The current constraint on spatial curva¬ 
ture from the Planck satellite, Vtk — 0.000±0.005 [54], implies that the curvature radius 
is at least 15idQ~^, and therefore that an observer will have causal access to regions of 
size C - Co < 1. 






• This procedure is then repeated for an ensemble of reference points xq, thereby tiling 
any constant field hypersurface with an ensemble of locally perturbed FRW patches. 

To compare our results with the CMB, we need to relate temperature anisotropies to 
the comoving curvature perturbation in the classical transition spacetime. From the vantage 
point of any observer, the perturbation associated with a classical transition covers the entire 
sky. In addition, since the perturbation is a pre-inflationary relic, we expect the affect to 
be primarily on the largest scales in the CMB. Different observers in the classical transition 
spacetime are sampled by considering different reference points xq. As we will see in the next 
section, for classical transitions there is a one-to-one mapping between xq and (^o- Below, we 
will use (^0 to label reference points. 

Working in the Sachs Wolfe approximation, the temperature fluctuation and TZ are 
related by [55] 

AT _ 7e(Xi,|xo) 

T 5 ’ ^ ^ 

where X\s is the intersection of the past light cone of an observer with the surface of last 
scattering in cartesian coordinates, a 2-sphere of radius 

^Is = 2y/Dfc(^o)^curv (3.16) 


where is the energy density in curvature today seen by an observer at (^o- In fho following, 
it will be convenient to compute the relative curvature in different patches. We therefore 
express relative to O.k{0) = ^k{^0 — 0): 

= (3.17) 

«o(^o) 


where = 0 is an arbitrary choice. The value of flk{0) is set by the details of the slow-roll 
inflation model to the future of the classical transition, and we will leave it arbitrary in what 
follows. 

Assuming that it is smooth (confirmed below), we can Taylor expand TZ about A = 0 
(equivalently ^ = ^o) and substitute with X = Rig cos 6 to obtain 


AT 

~Y~ 


1 

5 


2dlni^o\xo)nk{Co) cos2 e + 


(3.18) 


Expanding in spherical harmonics 

i m 


(3.19) 


we can identify 

020 = Y^^~^d‘^'^{^o\xo)Qki^o). (3.20) 

Substituting with Eq. 3.17, the contribution to the quadrupole of the angular power spec¬ 
trum ^ is therefore given by 


V2 

fifc( 0)2 


^ {dln{^o\xo)f 


Oo( 0 ) 

Oo(^o) 


^We report P2, defined as T>i = i{i H- l)C'^/ 27 r, where (airndi'm') — 


(3.21) 
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In summary, the leading observables in a classical transition spacetime are a correlated 
contribution to the spatial curvature (given by Eq. 3.17) and the quadrupole of the CMB 
temperature (given by Eq. 3.21), both expressed with respect to a model-dependent reference 
curvature ^2^(0). 

4 Simulation Results 

Let us begin by applying the pipeline described above to one of the AIP models. In Fig. 6 
we show a successful classical transition in an AIP model with the parameters A(/)i = 7.1 x 
10“^ Mp, /3i = 1.5, (j — 0.004Mp. The two expanding bubbles collide near x — l.OiL^^, 
creating a region of spacetime where the field rolls down the inflationary plateau. Choosing 
the constant field surface 00 = 0.15 M^, we show the mapping between a reference simulation 
co-ordinate xq and the reference anisotropic hyperbolic FRW coordinate <^o Fig. 7. The 
coordinate map is one-to-one with xq, anti-symmetric about the location of the collision, 
and diverges near the bubble wall. The wall reaches a constant comoving position in x, 
compressing the spatial extent of the slice in terms of (^o- 

The result for the observed temperature quadrupole and spatial curvature as a function 
of observer position is shown in Fig. 8. These plots show only half of the classical transition 
region, which is symmetric about the location of the collision at = 0. The anisotropy is 
maximized at the location of the collision at (^o = 0, falling to zero with increasing |(^o|- The 
magnitude of the anisotropy near = 0 is large unless the spatial curvature is relatively 
small. Comparing to the observed CMB quadrupole, the reference spatial curvature must be 
of order ^k{0) ^ 10“^ in order to produce a temperature quadrupole that is not too large 
(the observed CMB quadrupole is V 2 — 200 /iK^ [56]). The spatial curvature increases with 

plateauing to a constant value at large |(^o|- 

Because of the limited resolution and finite duration of the simulation, it is not possible 
to extend our description to arbitrarily large |(^o|- However, the quadrupole drops by roughly 
three orders of magnitude in the region we can compute. Extrapolating, we can expect that 
the universe becomes arbitrarily isotropic far from the centre of the classical transition, with 
fixed spatial curvature. If the curvature perturbation grew with distance from the centre, 
it would indicate that the surfaces of constant field were going timelike. Since the opposite 
occurs, there is every indication that the slices of constant field are everywhere spacelike, 
and therefore, that classical transitions produce spatially infinite universes. Most of the 
spatial volume on the constant field slices is far from the location of the collision, where the 
curvature perturbation vanishes. Therefore, most of the volume to the future of a classical 
transition is observationally indistinguishable from curved FRW. A similar result was found 
in Ref. [48], where it was found that in collisions between non-identical bubbles, infinite 
surfaces of constant field formed in the region to the future of the collision. The presence 
of domain walls appears to be the key ingredient in both examples, suggesting that inflating 
regions within domain walls are generically infinite in spatial extent. 

Turning to the BIP models, we first confirm that Normal, Oscillatory, and Marginally 
Repulsive geometries are produced by the simulation. Varying 02, which affects the depth 
of the intermediate vacuum as shown in Fig. 5, we indeed verify that each geometry can be 
produced in Fig. 9. This is a non-trivial check, as the classical transition produces a universe 
with a rolling scalar, and not a pure de Sitter space as was assumed in the thin-wall analysis 
of Ref. [52]. 
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Figure 6. A contour plot of (f) for a collision that produces a successful classical transition in the 
AIP model. The inflationary region of the potential is at (/> > .005 Mp. 
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Figure 7. Anisotropic hyperbolic coordinate, as a function of spatial simulation coordinate, Hpx. 
^ is anti-symmetric about the location of the collision {xHf = 1), and diverges near the bubble wall. 


In the Marginally Repulsive geometry, a lasting inflating region is produced. We use the 
pipeline described in the previous section to compute observables at different reference points; 
the result for a model with parameters A(p 2 = 3.5 x 10“^ ^2 — 9.865 and a — 0.004 

are shown in Fig. 10. As for the AIP model presented above, the temperature anisotropies 
are maximized near the centre of the classical transition region. However, the maximum 
quadrupole is 6 orders of magnitude larger than the AIP model in this particular example, 
requiring a far smaller reference curvature to be consistent with the observed quadrupole 
if we were to inhabit this region. In contrast to the AIP model, the spatial curvature has 
no signs of plateauing in the region that the simulations can accurately probe. The fall-off 
of the quadrupole suggests that there is no obstruction to having spatially infinite hyper 
surfaces. However, the growth of curvature suggests that the Universe near the location of 
the collision will appear quite different from regions far from the location of the collision. 
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^0 ^0 

Figure 8. The temperature quadrupole P 2 in (left) and spatial curvature (right) observed 
at different reference points ^0 for the AIP model with parameters: A0i = 7.1 x 10-4 M^, A = 1.5, 
cr = 0.004 Mp and Ax = 1.9i7^^. From the map shown in Eq. 7, ^0 = 0 can be identified with the 
location of the collision. These plots depict only half of the classical transition spacetime. The other 
half is obtained by mirroring the plots across the right vertical axis. The anisotropy is largest near the 
centre of the classical transition, falling off with distance. The spatial curvature appears to asymptote 
to a constant value a few times larger than its value in the centre of the classical transition, indicating 
that fewer e-folds of inflation occur at the edges of the classical transition than in the centre. 
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Figure 9. Bubble collision simulations for various (32. For small P 2 a Normal geometry is produced 
(left panel), for intermediate P 2 a Marginally Repulsive geometry is produced (centre panel), and for 
a large /32 an Oscillatory geometry is produced. See Fig. 1 for a detailed summary of the various 
possible classical transition geometries. 


The dramatic growth of curvature could imply negligible inflation far from the location of the 
collision, although we cannot determine the extent of this effect given the finite resolution of 
our simulations. 

4.1 Scanning the parameter space of the AIP model 

As described in Sec. 3.1, for the AIP model we explore slices through parameter space which 
preserve the radius R and amplitude 0amp of the bubble walls. In addition, we sample different 
bubble separations Ax. Qualitatively, the results over the range of parameters sampled are 
identical to those described in the previous section. Specifically, the temperature anisotropy 
peaks at the centre of the classical transition, falling to zero with increasing j^ol, and the 
spatial curvature increases to a constant at large |^o|- Variations in the initial separation Ax 
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Figure 10. The temperature quadrupole P 2 in (left) and spatial curvature ftk (right) observed 
at different reference points ^0 for a marginally repulsive geometry in the BIP model with A02 = 
3.5 X 10“"^ Mp, /32 = 9.865 and a = 0.004 Mp. The anisotropy is largest near the centre of the classical 
transition, falling off with distance. The spatial curvature increases with |^o|: indicating that inflation 
is significantly disrupted far from the centre of the classical transition. 


and parameter combinations that keep the amplitude 0amp fixed but vary the initial radius R 
control the centre of mass energy in the collision (in flat space, the lorentz factor is given by 
Eq. 2.1). Parameter combinations that keep the radius R and initial separation Ax fixed but 
vary the amplitude (/)amp are independent of the centre of mass energy, and directly related 
to the underlying scalar field lagrangian. For a specific potential and initial separation Ax, 
the value of V 2 and are correlated at each reference point. We can therefore represent 
the predictions of a specific model as a line through the V2-^k plane. 

In Fig. 11 we show the T>2-^k plane for models where we vary Ax and i?, which control 
the centre of mass energy of the collision. In the left panel, we see that increasing Ax 
generally leads to larger V 2 and Vtk- The fall-off of V 2 and plateauing of Vtk at large |^o| can 
be seen as one goes from left to right along each curve. A particular microphysical model 
(e.g. scalar field potential) yields an ensemble of predictions given by the set of curves in the 
V2-^k plane for all Ax. For a specific scalar field lagrangian, bubble collisions will sample 
all possible values of Ax with probability Pr oc sin^ Ax [57] . In the centre and right panel 
of Fig. 11, we show the predictions in the V2-^k plane for varying initial radius i?, keeping 
both Ax and 0amp fixed. In general, smaller initial radii (and therefore a higher centre of 
mass energy) yield larger V 2 and Vtk- Again, the fall-off of V 2 and plateauing of Vtk at large 
|(^o| can be seen as one goes from left to right along each curve. In contrast to the initial 
separation, the initial radius is a deterministic parameter fixed by the scalar field potential. 
In conclusion, the observational signatures of a classical transition clearly depend on the 
centre of mass energy of the collision. More specifically, collisions with a higher centre of 
mass energy lead to larger predicted values for V 2 and Q^k- 

In Fig. 12, we show the V2-^k plane for models where the initial separation Ax and 
radius R are held fixed, but (/)amp varies. Due to computational limitations, we are only 
able to sample a relatively limited range of 0amp- In the left panel, we vary and A(j)i 
while holding a fixed at cr = 0.006 M^. Here, there is an approximately linear decrease in 
the maximum value of V 2 and an increase in the maximum value of Q^k- The increase in 
the maximum value of Q^k is consistent with the free-passage approximation, since the field 
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Figure 11. Correlation between V 2 in iiK‘^ and iii the AIP model. In the left panel, we show 
variations in the initial separation Ax. Curves correspond to Ax = 1.1 (blue), and Ax = 1.9 (green). 
Ax = 2.3 (red). In the centre panel, we show variations in initial radius R from the Acjyi-a sector 
of the potential. Curves correspond to HpR = .246 (blue), HpR = .477 (green), and HpR = .728 
(red). In the right panel, we show variations in initial radius R from the /3i-a sector of the potential. 
Curves correspond to HpR = .22 (blue), HpR = .53 (green), and HpR = .79 (red). 
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Figure 12. Correlation between V 2 in iiK‘^ and fhe AIP model. In the left panel, we show 

variations in the instanton amplitude 0amp from the ^i-Acjyi sector. Curves correspond to (/)amp = 
2.3 X 10“^ Mp (blue), 0amp = 2.8 x 10“^ Mp (green), and 0amp = 3.3 x 10“^ Mp (red). In the centre 
panel, we show variations in the instanton amplitude 0amp from the A(/)i-cr sector. Curves correspond 
to 0amp = 1.0 X 10“^ Mp (blue), 0amp = 2.2 X 10“^ Mp (green), and 0amp = 3.6 x 10“^ Mp (red). In 
the right panel, we show variations in the instanton amplitude 0amp from the I3i-a sector. Curves 
correspond to 0amp = 0.9 x 10“^ Mp (blue), 0amp = 1.5 x 10“^ Mp (green), and (/)amp = 2.5 x 10“^ Mp 
(red). 


will enter inflation later for a larger 0amp- The decrease in the maximum value of P 2 can be 
understood as due to the delayed horizon crossing of 7^, which leads to additional decay. In 
the A(/)i-cr and /3i-cr sectors, there is no apparent trend associated with varying 0amp- This 
is due to the fact that as a is varied, the field value at which slow-roll begins is varied. 

Returning to the form of the potential (see Fig. 3), we see that if a is large enough then 
a period of field evolution will happen before the collision. It was found in Ref. [17] that field 
evolution, in addition to the barrier width, can contribute to the field excursion produced 
by a bubble collision. To explore this behaviour in more detail, we perform simulations in 
which we vary a. For example, choosing A(/)i = 0.0011, A (/)2 = 0.0007, — 2.0 and ^2 — 3.0 

a classical transition is obtained for a < .016. In this example, the distance the field is 
displaced by the collision is five times iarger than the ampiitude of the instanton, which is 
^ .003. 

Adjusting a such that we approach the boundary where ciassicai transitions are aiiowed, 
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Figure 13. Growth of lumpiness as the second bump approaches classical transition boundary. The 
potential parameters are set to A(/)i = 0.001061 Mp, A (/)2 = 7.1 x 10“^ Mp, Pi = 2.0 and P 2 = 3.0. In 
the left panel, the shift is set to a = .007, where there are few wall modes excited. In the right panel, 
the shift is set to a = .016, just at the edge of the range where classical transitions are possible. The 
wall modes are visible as oscillations in the comoving curvature perturbation. 


the comoving curvature perturbation develops additional structure. In Fig. 13 we show 
which develops an oscillatory structure as the classical transition boundary is approached. 
These features are produced by the excitation of the internal degrees of freedom of the post¬ 
collision domain walls. Such wall modes were also observed in Ref. [48] for the collision 
between non-identical bubbles. Wall modes add extra structure to the comoving curvature 
perturbation beyond the quadrupole. We leave an exploration of the phenomenological con¬ 
sequences of wall modes to future work. 

4.2 Exploring the BIP model 

A marginally repulsive geometry is produced for a specific relation between the vacuum 
energies and the energy density at the onset of slow-roll inflation, as described in Sec. 2. 
The vacuum energies of the two minima are unambiguous. However, we must impose some 
prescription for measuring the energy density at the onset of slow-roll. We define the onset of 
slow-roll as the field value (psp along the first space like surface of constant field; an example 
is shown in Fig. 14. 

In Fig. 15 we show the V2-^k plane in the two sectors (cr — P 2 and A (/)2 — P 2 ) of 
parameter space that give marginally repulsive classical transition geometries. Comparing 
with the results from the AIP model shown in Figs. 11 and 12, there are a number of notable 
differences. First, the size of both V 2 and is orders of magnitude larger in the BIP 
models. In addition, there is no clear asymptote of at large ^o- Looking at the trends 
associated with the field value (j)^^ at the onset of slow-roll, we see that the signature grows 
with increasing This makes intuitive sense, as more violent initial perturbations will 
delay the onset of inflation. Finally, we have explicitly verified that the energy density at 
the onset of inflation is above the energy density at the intermediate potential medium, but 
obeys the bound Eq. 2.2. 

5 Conclusion 

In this work, we investigated the possibility of observationally determining whether our ob¬ 
servable Universe arose from a classical transition produced by the collision of two bubbles. 
We consider single-field models with a potential that has three local minima, as depicted in 
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Figure 14. An example of countour plot for the second model in the ^2 — o' sector. The magenta line 
indicates the first space like constant-field hyper surfaces, which we take as the beginning of inflation. 






Figure 15. Correlation between V 2 in fiK^ and Clk for the BIP model. In the left panel we sample 
different values of 0sp by varying parameters in the (7 — ^2 sector. Curves correspond to 0sp = 0-01 
(blue), 0sp = 0.018 Mp (green), and (/)sp = 0.025 Mp (red). In the right panel we sample different values 
of (/)sp by varying parameters in the A02 — P 2 sector. Curves correspond to 0sp = 0.0130 Mp (blue), 
0sp = 0.0121 Mp (green), and (/)sp = 0.0114Mp (red). 


Fig. 1. To have a viable model of our observable Universe, we additionally require a slow-roll 
plateau leading to the lowest energy minimum. There are two classes of models that could 
give rise to our observable Universe. The first class of models (AIP models) produce Normal 
classical transition geometries. These models have three minima with sequentially decreasing 
energy density, and the collision produces an ever-expanding region containing the lowest en¬ 
ergy state. The second class of models (BIP models) produce Marginally Repulsive classical 
transition geometries. In the BIP models the classical transition produces an inflationary 
epoch of a higher energy than that of the bubble interiors (but obeying the bound Eq. 2.2). 

In both models, we have found that classical transitions produce approximately FRW 
Universes with negative spatial curvature that are (extrapolating our simulation) infinite 
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in extent. The classical transition geometry is inhomogeneous, but retains approximate 
planar symmetry. We expand the classical transition geometry into local approximately FRW 
patches characterized by the local spatial curvature and the comoving curvature perturbation. 
The leading observables in such a Universe are a CMB quadrupole and negative spatial 
curvature. 

Our main results are presented in Figs. 11, 12, and 15 which show the the V2-^k plane 
for varying parameters in the models under consideration. All of our results are presented as 
a function of an arbitrary reference curvature ^2)^(0), which depends on the details of slow-roll 
and reheating. In the AIP model, the spatial curvature increases to a constant with distance 
from the centre of the classical transition. The locally observed quadrupole drops to zero with 
distance from the centre. The locally observed quadrupole also drops to zero with distance 
from the centre of the classical transition in the BIP model, however, the spatial curvature 
appears to increase without bound. 

The magnitude of the locally observed quadrupole and negative spatial curvature de¬ 
pends both on the centre of mass energy of the collision (fixed by the initial separation Ax 
and initial radius R of the colliding bubbles) and the scalar field lagrangian. In both the AIP 
and BIP models, decreasing initial radius or increasing initial separation generally leads to 
a larger locally observed quadrupole and curvature. Varying parameters in the scalar field 
lagrangian, we choose parameter combinations in the AIP model that fix the initial radius of 
the bubble, but change the total field excursion (/)amp between the instanton endpoints. We 
find that an increasing amplitude leads to a larger locally observed quadrupole and curvature. 
We have also found that motion of the field after tunnelling can contribute significantly to 
the distance in field space that the classical transition can reach. Finally, we have found that 
the amplitude of the locally observed quadrupole and curvature is generally far larger in the 
BIP models than the AIP models. 

Although we have not investigated how changes in the inflationary potential will affect 
our results, in analogy with previous results on bubble collisions [17], we can expect that 
small-field models of inflation to be more susceptible to disruption than large-held models 
of inhation In the moments after a collision, the details of the potential are irrelevant. 
Therefore, we can compare the held excursion predicted by the free passage approximation to 
the size of the region in held space over which inhation occurs. For a hxed barrier structure, 
the held excursion produced by the classical transition is hxed, and so small-held models will 
experience a larger fractional displacement than large-held models. The held excursion during 
inhation is related to the tensor to scalar ratio r by the Lyth bound [58]: A(l)lMp ^ y^r/0.01. 
For a hxed barrier structure, we therefore expect that the duration of inhation will be affected 
more severely in models with small r. The duration of inhation in a given model (modulo 
the details of reheating) determines the reference curvature riyt(O), which in turn determines 
the magnitude of observables. Therefore, inhationary models with small r can be expected 
to yield a larger spatial curvature and a larger locally observed quadrupole. 

Will we ever be able to test whether or not our observable Universe was created by a 
classical transition? The fact that the observables are only curvature and a quadrupole makes 
it difficult to conhrm whether or not a classical transition in a single held model occurred in 
our past. Curvature can be produced in many scenarios, and the temperature quadrupole has 
an enormous error bar due to cosmic variance. Ruling out models with a classical transition 
requires an observation of spatial curvature. Current observation constrains curvature at 

^ Small field models are defined as those in which inflation takes place over a sub-planckian field range, 
while large field models are defined as those in which inflation takes place over a super-planckian field range. 
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the level of Qk — 0.000 ± .005. The confusion limit with super horizon perturbations is 
^ ±10“^, which represents a lower bound on what level of primordial curvature could 
ever be observed. An observation of positive spatial curvature at a level greater than this 
would rule out a classical transition in our past. If an observation of negative spatial curvature 
is made, one can look at the correlated quadrupole in a given model, and compare to the 
observed value. However, in the AIP models, a particular value of the spatial curvature can 
be correlated with arbitrarily small values of V 2 due to the fact that curvature asymptotes 
to a constant at large distance from the centre of the classical transition, while the locally 
observed quadrupole drops to zero. In the BIP models, the curvature grows with distance 
from the centre while the quadrupole shrinks. Therefore, for a fixed ^k{0) it is conceivable 
that one could rule out some portion of model space based on the non-observation of a 
large curvature or a large quadrupole, since one or the other will be predicted depending on 
the location of the observer. However, ^k{0) is dependent on the number of e-folds in the 
region that underwent a classical transition, which is completely unconstrained. This makes 
it impossible to conclusively rule out a classical transition in our past. 

Signatures in multi-field models could be more varied than single field models, for ex¬ 
ample giving rise to isocurvature perturbations or leading to new sources of instability (see 
e.g. [59-61]). In addition, a deeper exploration of the role of wall modes is warranted, and 
may yield a more distinctive signature than a CMB quadrupole. This may make it easier 
to constrain models where our Universe arose from a classical transition. In future work, 
we plan to explore the phenomenology of bubble collisions and classical transitions in multi 
field models. It would also be desirable to embed a model of classical transitions in a more 
realistic theory, for example string theory. Indeed, classical transitions may be expected 
in the context of flux compactifications [35]. Regardless, this paper illustrates that models 
of eternal inflation can make predictions for real observables, adding further credence to a 
scientific theory of the multiverse. 
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